Extracting chemical energy by growing disorder: Efficiency at maximum power 
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We consider the efficiency of chemical energy extraction from the environment by the growth 
of a copolymer made of two constituent units in the entropy-driven regime. We show that the 
thermodynamic nonlinearity associated with the information processing aspect is responsible for 
a branching of the system properties such as power, speed of growth, entropy production, and 
efficiency, with varying affinity. The standard linear thermodynamics argument which predicts an 
efficiency of 1/2 at maximum power is inappropriate because the regime of maximum power is 
located either outside of the linear regime or on a separate bifurcated branch, and because the usual 
thermodynamic force is not the natural variable for this optimization. 

PACS numbers: 05.70.Ln,05.70.-a,05.40.-a 



I. INTRODUCTION 

Carnot efficiency is one of the cornerstones of thermo- 
dynamics since it leads to the definition of entropy and 
the Second Law of thermodynamics. It expresses a fun- 
damental limitation on how much work can be extracted 
from a heat flow. A less studied but arguably more rele- 
vant question for many isothermal chemical and biologi- 
cal processes is how much chemical energy a system can 
extract from its environment by increasing the system's 
configurational entropy. Thermodynamics does, in fact, 
also prescribe a limit, even though at first sight it ap- 
pears to be almost trivial: the energy extracted by such 
an isothermal transfer can be carried out with 100% effi- 
ciency. However, there is a crucial additional condition, 
namely, that this efficiency can only be reached - just 
as in the case of Carnot efficiency - by a reversible, in- 
finitely slow process. Hence 100% efficiency is achieved 
for a process with zero power output. The question of 
efficiencies at finite power should thus be addressed. 

In the context of thermal machines, a straightfor- 
ward analysis based on linear irreversible thermodynam- 
ics teaches us that, as one moves away from the reversible 
regime, the power goes through a unique maximum, and 
that the efficiency at this maximum is, at most, 50% 
0, The same argument can easily be extended to 
the transformation between different forms of chemical 
energy. However, the above prediction may not apply 
for several reasons. First, the point of maximum power 
does not necessarily lie in the linear regime. Second, 
thermodynamic nonlinear effects can give rise to bifur- 
cated branches. Finally, the above-mentioned efficiency 
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is attained upon maximization with respect to the ther- 
modynamic force associated with the power generating 
flux. While this is a natural set-up in many problems, it 
may not always correspond to the relevant scenario. 



In this paper, we investigate the efficiency of a chem- 
ical entropy-driven process of capital importance in bio- 
physics, namely, copolymer synthesis [!, 0, H, H[ , see also 
0,1a]. As exemplified by the copolymer DNA, guardian 
of genetic information, such processes are essential for 
biological information processing. We will show that 
the above-mentioned complications are present in this 
generic model. In particular, the thermodynamic nonlin- 
earity associated with the information processing aspect 
is responsible for a branching of the systems properties 
such as power, speed of growth, efficiency and entropy 
production, as one varies the affinity. Furthermore, the 
regime of maximum power is located either outside of the 
linear regime or on the separate bifurcated branch. Fi- 
nally, it turns out that the thermodynamic force is not 
a natural control variable in the present model. While 
our (exact) analysis is carried out for the simplest pos- 
sible model, namely, copolymer synthesis with two con- 
stituent building blocks, our findings suggest that chem- 
ical information processing usually operates in the far- 
from-equilibrium regime, with unique features due to the 
entropic contribution. 



In Sec. |TT]we present the basic thermodynamic formu- 
las that define our system. In Sec. IIIII we present the 
detailed kinetic description of our model, whose results 
are discussed in detail in Sec. IIVI In particular, it is here 
that we exhibit the correct and unexpected results for 
the efficiency at maximum power, results that arise en- 
tirely from the nonlinear nature of the problem. A brief 
recapitulation is presented in Sec. [Vj 
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II. THERMODYNAMICS 



thus be written in the familiar bilinear form 



We begin with some well-known relations for isother- 
mal systems. Consider a spontaneous chemical process 
involving particles of different types labeled by j, with 
corresponding particle number Nj and chemical potential 
fij . The system is in thermal and mechanical equilibrium 
at temperature T and pressure P. The total Gibbs free 
energy 

G = U + PV - TS = H - TS = IMNi (1) 

i 

evolves toward a minimum value, so that dG < 0. Alter- 
natively, to characterize the evolution of the isothermal 
system we write 

dS = diS + d e S 

TdiS = — HjdiNj 
j 

Td e S = dH (ijdeNj . (2) 



We have separated the total entropy change into two con- 
tributions. The first one, diS, is the always-positive part 
of the entropy change, called the internal entropy produc- 
tion. The other is the contribution to the entropy change 
due to exchange processes between the system and its en- 
vironment, and can be positive or negative. Associated 
with these contributions, we have written the change in 
he number of particles of type j as 



dNj = d e Nj + diNj, 



(3) 



where the first contribution is due to exchange of parti- 
cles with the environment, and the second is the internal 
change caused by the chemical reaction. We take the sys- 
tem to be closed, i.e., it exchanges only energy but not 
particles with the environment, so that d e N = 0. These 
definitions lead to consistency between the statements 
that the system evolves toward a minimum in the Gibbs 
free energy and that the internal entropy production of 
this chemical system has to be positive [9j, that is, 



diS = dS 



dH 
~T~ 



Ei NdNj dG 



T 



T 



> 0. (4) 



Obviously, for a reversible transformation with zero in- 
ternal entropy production diS = 0. 

We now turn to the simplest scenario of copolymer 
synthesis. The system consists of a bulk-phase containing 
two types of monomer units, f and 2, which can attach 
or detach at the endpoint of a single linear copolymer. 
We identify the four states 1/, fc, 2f and 2c. Here jf 
represent free bulk monomers and jc represent monomers 
attached to the copolymer. Since the number of each 
type of monomer is conserved, one has dN% / = —dN\ c 
and dN2f = —dN 2c . The entropy production (0J can 



djS _ 



, dN lc 



(Vif - Mic)-^- + (M2/ - M2c) 
-A 2 J 2 , 



dN 2l 
dt 



(5) 



with the affinities Aj = (/x,-/ 
fluxes Jj = dNjc/dt. 

In view of the relation diS 
the entropy production as 



Si 



«2 



T 



[ijc) and the conjugate 
dS — dH/T, we rewrite 

)J 2 +D{J l + J 2 ). (6) 



Here hj is the change of enthalpy per monomer upon 
transfer from the bulk to the copolymer. The crucial 
point, which has been discussed in detail in the litera- 
ture [1, H, H 0] , is to realize that the average change of 
entropy upon transfer of a monomer from the bulk to 
the copolymer contains two contributions. One is the 
monomer entropy, Sj , due to the change in the monomer 
degrees of freedom and in the monomer internal structure 
between the free monomer in solution and the monomer 
inside the copolymer. The other is the configurational 
entropy denoted by D, due to the change in the informa- 
tion stored in the polymer sequence that occurrs when a 
monomer is added to the copolymer. It is given by the 
Shannon entropy 



D = - lim - V P^lnP^ 

l^oa I ^— ' 



(7) 



where / is the copolymer length in monomer units and P u 
is the probability of a copolymer with monomer sequence 
lo. In the absence of correlations, the Shannon entropy 
is expressed solely in terms of the monomer abundance 
probabilities pi — p and p 2 — 1 — p, 



D = —p\np — (1 — p) ln(l — p). 



(8) 



For simplicity, we further assume that monomer entropy 
and enthalpy changes upon transfer of a monomer from 
the bulk to the copolymer and vice versa have the same 
value for both monomers, that is, 



hi _ h# 

T ~ si - T 



s 2 . 



(9) 



We henceforth call Te the monomer "free enthalpy." In- 
troducing the net speed of growth of the copolymer, 
v = Ji + J 2l the entropy production can finally be written 
as 



Si = Av > 0, 
where the total affinity is given by 
A = D-e. 



(10) 



(11) 



The expression (|I0[) for the entropy production in the 
steady state regime of the growing copolymer has been 



3 



derived in Refs. 0, 0, S Q ■ It is interesting to realize that 
the affinity is not an obvious control parameter due to its 
dependence on D which is in turn a nontrivial function of 
e. Only e can be easily controlled externally by changing 
the concentration of the monomers in solution. 

The power at which e, the free enthalpy divided 
by temperature, is extracted from the surroundings by 
copolymer growth is given by 



ev = (D - A)v, 



(12) 



The efficiency rj of the process is defined as the ratio of 
this power over the cost Dv of the entropy growth per 
unit time, 



V 



ev 
~Dv 



D-A 
D ' 



(13) 



In the reversible limit with A,v — > 0, the efficiency of 
the process becomes optimal, rj = 1, but the extracted 
power goes to zero. The standard prediction from linear 
thermodynamics that arrives at an efficiency of 50% at 
maximum power is obtained upon expanding the velocity 
in terms of the affinity, v = LA, with L the linear re- 
sponse coefficient. Within this approximation the power 
becomes P » L(D — A) A. Note that this power attains 
its maximum for A = D/2 with the corresponding effi- 
ciency T] = 1/2, if we assume that D is kept constant. 
However, below we will investigate the more natural op- 
timization with respect of e, since this is the natural and 
easily controllable variable related to the free enthalpy 
flux. Whatever control variable is used, we will see in 
Sec. IIVI that the true maximum is beyond the reach of 
this linear expansion (and even of a nonlinear continua- 
tion of this expansion) . 



this occurs with probability p. The net rate of attach- 
ment for 2 is similarly given by fc +2 — k-%{l — p)- We 
thus conclude that 



fe+i - k-ip 



P 



k +2 -k- 2 (l-p) l-p 



(15) 



The solution of the resulting quadratic equation for p 
reads 



P 



y/a 2 - 4(fc_i - fc_ 2 )fc^ 
2(A;_ 1 -fc_ 2 ) 



(16) 



where a = k+\ + fc+ 2 + k— l — k—2- By a similar argument 
we find that the speed of growth of the copolymer, given 
by the rate of attachment fc+i + minus the rate of 
detachment fc_ip + fc_ 2 (l — p), is given by 



v = k +1 - k_ x p + k +2 - fe_ 2 (l - p). 



(17) 



We note from Eqs. ©, (10]), (TTJ, ([IB]) and fI7|) that 
equilibrium, v — and A — 0, occurs at e = In 2 with 
p = 1/2 and D = In 2. For smaller (larger) values of e, 
A > (A < 0) and the copolymer is synthesized (de- 
graded), i.e., v > (v < 0). Of specific interest to us is 
the surprising regime of entropy-driven growth, A > 
and v > 0, but with e > @, 0, H, @. Under the 
simplifications assumed in our model, this occurs when 
< e < In 2. Monomers are pumped uphill against the 
free enthalpy barrier e > under the influence of the 
entropic contribution D to the affinity. The power P (en- 
thalpy per unit time) extracted from the copolymeriza- 
tion dynamics is positive in this entropy driven regime, cf. 
Eq. (|12p , with corresponding efficiency given in Eq. (11 



III. KINETICS 

We now turn to the detailed kinetic description of the 
copolymerization process, which will allow us to identify 
the expressions for v and p in the context of a full nonlin- 
ear analysis. Let us call k + j and the rates of insertion 
and removal, respectively, of monomer j = 1,2. Because 
the free enthalpy of the monomers has been assumed to 
be the same, the ratios of the reaction rates are given by 



The fraction p of monomers of type 1 present in the 
copolymer in the regime of steady growth can be deter- 
mined by the following self-consistency argument. The 
ratio p/(l — p) of the number of 1 versus 2 monomers in 
an ensemble of copolymers has to be equal to the ratio 
of their net rates of attachment. For monomer 1, this 
net rate is the pure rate of attachment, fc+i, minus the 
rate of detachment, which is — k—ip. The factor p arises 
from the fact that detachment is only possible when the 
monomer at the tip of the copolymer is of type 1, and 



IV. EFFICIENCY AT MAXIMUM POWER 

To study the main question of interest, namely, the 
regime of maximum power and its corresponding effi- 
ciency, we choose convenient variables. We note that 
the model is described by four kinetic constants, but the 
latter are not independent since they obey the relation 
Eq. (|14[) . Furthermore, one of them can be set equal to 1 
by an appropriate choice of the time unit, e.g., k + \ = 1. 
As the remaining two degrees of freedom, we choose e 
and k + 2- We then have explicit functional expressions 
for all the other quantities = e e , fc_ 2 = e e fc +2 , p 
and v, cf. Eqs. (jTHJ) and (|17p. and hence also D, A, P, 
Si, and r), see Eqs. (HI]), Q2J|, O, © and OH- Other 
relations between, for example, P and A, can then be ob- 
tained by parametric elimination. The quantities P, v, 
■q, A, and Si can easily be calculated numerically. The 
results are summarized in Figs. [T][2][3l and|4] We next 
turn to a discussion of these figures, supplemented with 
corresponding analytic calculations. 

The information contained in Figfl] is detailed in the 
caption. The most striking feature in this figure is the 
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FIG. 1: (Color online) The full thick curves represent the 
parametric dependence on < e < In 2 of the power P, the 
copolymerization speed v, and the efficiency 77 on the ordinate 
and the affinity A as abscissa. The point e = In 2 is located 
at the origin of the axes for P and v and at r\ = 1 and A — 

for 77. The small dots along the curves are separated by 
Ae = (In 2)/14 to indicate how fast e changes along the curves. 
Different thick curves correspond to different choices of k+2, 
with fc+i < k+2 < oo. Without loss of generality we set fc+i 

1 (time rescaling). The thin dashed curves intersect the thick 
curves where the value of e corresponds to maximum power 
with respect to e. The curves in the inset in the P plot and 
in the upper inset of the v plot correspond to k+2 = 65, 85, 
105, 121, 160, and 200. The dashed curves in the lower inset 
of this plot represent the linear response predictions v — LA 
for k +2 = 1,2,11. 
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FIG. 2: (Color online) Efficiency 77 and power P corresponding 
to the maximum power denoted simply by e in the figure for 
different values of k+2- We have set fc+i = 1. 
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FIG. 3: (Color online) Derivative of the affinity with respect 
to e and affinity (in the inset) as a function of e. The different 
curves correspond to k+2 ~ 0, 11, 65, 85, 105, 121, 160, 200, 
1001, and 00 and fc+i = 1. 



existence of two different branches for the power and ve- 
locity in terms of the affinity. The transition between the 
two branches occurs when 



v) _ d(P, v) (dA 
dA de \ de 



(18) 



diverges. (P, v) indicates P or v. Since d(¥,v)/de is 
an analytic function of e, the new branch appears when 
dA/de touches zero. As long as the latter quantity re- 
mains positive, which is the case for k + 2 smaller than 
a certain critical value, cf. Fig. [3l the power and ve- 
locity can be seen as a true function of A. Branch- 
ing takes place at the critical point, characterized by 
dA/de = d 2 A/de 2 = 0, resulting in fc +2 « 84.33 and 
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e = 0.088, see again Fig. [3] For values of fc +2 larger than 
this critical value, power is no longer a proper function 
of A, as two branches appear, with two values of (P, v) 
for a given value of A. While along the linear branch and 
its continuation the affinity decreases with e, the affinity 
increases with e on the new lower branch, cf. the inset 
in Fig. [3] This remarkable result implies that we can 
approach low values of affinities via a nonlinear branch 
which is distinct from the branch predicted by linear re- 
sponse theory and its continuation. We note that the en- 
tropy production itself becomes a bi-valued function in 
terms of the affinity, as can be seen in Fig. 0J Naively, one 
would expect that entropy production and affinity both 
provide consistent measures for the distance from equi- 
librium. This is clearly not the case in the present model, 
where the entropy production is a decreasing function of 
the affinity on the upper nonlinear branch. In particular, 
for very large values of fc+2 one finds that the entropy 
production becomes very large while the affinity goes to 
zero. We conclude that the affinity is not a reliable mea- 
sure for the distance from equilibrium. 



The efficiency thus becomes 
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FIG. 4: (Color online) Same type of plot as Fig. [1] but for 
entropy production. 

To explore the region close to equilibrium and, in par- 
ticular, the linear response regime, we write e = In 2 — Se 
and expand in powers of 5e. From Eq. ([5]) with Eqs. UJ 
and (11411 we find 



where 



D = ln2- a Se 2 + 0(Se 3 ), 



(fc+i - k +2 f 



2(fc +1 +fc +2 ) 2 ' 
For the affinity, we find from Eq. (fTTj) that 
A = 6e-a8e 2 + 0(5e 3 ). 



(19) 



(20) 



(21) 



iS + InV^ 3 ) 
A 



= 1 -hT2 + ^ 



(22) 



This linear dependence of the efficiency on the affinity 
close to equilibrium is clearly identified in the upper 
left region of the third affinity plot in Fig. [TJ while the 
corresponding behavior of the affinity in terms of e, cf. 
Eq. (|21[) . is observed in the lower left region of the inset 
of Fig. [3] In this regime close to equilibrium we find the 
standard linear response relations 



= LA + 0{Se 2 ) 

= In 2 LA + 0{6e 2 ), 



(23) 
(24) 



with the Onsager coefficient given by L = 
4/c+ifc+2/(fc + i + k +2 ), cf. the lower left regions of 
the power and speed plots in Fig. [TJ Note also that the 
Onsager coefficient becomes independent of fc + 2 in the 
limit k + i <C fc+2, where L = 4fc+i. 

We have seen that linear response predicts an efficiency 
at maximum power of 50%. However, as announced ear- 
lier, this result is not correct. This is seen in Fig. [1] or in 
Fig. [21 where the affinity is clearly above the value 1/2 in 
the regime "closest" to equilibrium. The explanation is 
that maximum power occurs beyond the reach of linear 
response theory, as can clearly be seen in the lower inset 
of the v plot in Fig. [TJ where the linear response curves 
(dashed lines) become inaccurate at maximum power. 
Furthermore, we note that the point of maximum power 
moves onto the nonlinear branch as fc+2 grows, now oc- 
curring at decreasing values of A. So, even though we are 
approaching a regime of low power output with decreas- 
ing affinity, we do so via the nonlinear branch, where the 
prediction of linear response theory utterly fails. The 
main conclusion is that, while there is indeed a regime 
of linear response, it is unable to describe the region of 
maximum power, which always occurs outside the regime 
of validity of the linear law. 

To complete our analysis, we explore in detail the lim- 
iting cases k+2 — > fc+i and fc+2/fc+i — > oo. For trans- 
parency, we explicitly retain fc+i instead of setting it 
equal to unity. In the limit where k+2 fc+i, we find 
that 



1 



v = fc+i(2-e e ), L> = ln2. 



(25) 



This leads to an efficiency r\ = e/ln2 = 1 — Aj In 2, as 
observed in Fig. [TJ In this limit, the value of e leading 
to maximum power is obtained as the solution of the 
transcendental equation 2e~ e — e = 1, namely, e « 0.375. 
At maximum power we thus get 



0.204fc+i and n ps 0.541, 



(26) 



as seen in Figs. [TJ and [5J As an immediate consequence, 
we also find v w 0.545fc+i and A ss 0.318, as observed in 
Fig.rT] 
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In the limit k+ 2 — * oo, where 



p=l 



k+ie' 



(2e" 



1) 



(e- - 1) ' 
L» = c- £ e-(l-e- e )ln(l-e- e ), 



(27) 



the efficiency reads 77 = — e/[c~ e e + (1 — c~ e )ln(l — 
e~ c )]. The numerical results of Fig. Q] suggest that max- 
imum power in this limit occurs for e very close to 
zero. We therefore expand the velocity around e = 
and find v = -y/k+ik +2 + (k+i + k +2 )e/2 + 0(e 2 ). 
Using Eq. (|12|) . we find that maximum power occurs 
at e = — -\/fc+ifc + 2/(fc + i + fc+2), resulting in P = 
yJk^ik+2 1 '{k+i + k +2 ). For k +2 — » 00 the latter becomes 



1, 



(28) 



as observed in Figs. [T]and[21 Similarly, by expanding rj 
to first order around e = and using the value we found 
for e at maximum power, we find that 



»7-»0, 

as observed in Figs. [T] and [2] 



(29) 



V. CONCLUSIONS 

Using a simple model of copolymerization, we have 
shown that free enthalpy can be extracted from the en- 
vironment in response to the entropic force correspond- 
ing to the information stored in a growing copolymer 



sequence. The thermodynamic nonlinearity associated 
with the information processing aspect is responsible for 
a branching of the dependence on the affinity of system 
properties such as power, speed of growth and efficiency. 
The nonlinear regime occuring after the branching is par- 
ticularly surprising since the entropy production keeps in- 
creasing even as the affinity begins to decrease. We iden- 
tified a regime of linear response where the efficiency of 
the energy extraction is optimal (equal to 1), but where, 
as usual, the power output goes to zero. Considering in- 
stead the efficiency at maximum power, we found that 
the universal prediction of linear response theory (effi- 
ciency equal to 1/2) is inappropriate for this model. The 
reason is that the copolymerization generating maximum 
power occurs far from equilibrium in a region not accessi- 
ble to linear response theory. Our results suggest a possi- 
ble self-powering mechanism for nonequilibrium systems 
that can extract chemical energy from their surroundings 
by growing their internal structural information. 
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